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Discussion of 
PAVLOVSKY’S THEORY FOR PHREATIC 
LINE AND SLOPE STABILITY 


by K.P. Karpoff, A.M., ASCE 
(Proc. Sep. 386) 


HAMILTON GRAY,! M. ASCE.—This discussion is limited to comments on 
the second portion of Mr. Karpoff’s paper which appears under the heading, 
“Stability of Slopes”. Certain modifications appear to be appropriate in the 
interest of accuracy of analysis, although the actual numerical discrepancies 
which are involved do not appear to be serious. The deterivration of the face 
of a slope by seepage water in the manner contemplated by Mr. Karpoff occurs 
not only in artificial earth dams and embankments, but also in natural slopes 
where the face of the slope intercepts the natural position of the ground water 
surface. Usually this process of deterioration is observed after a slope has 
been newly formed or altered in connection with construction, since the deter- 
ioration, if allowed to continue unchecked, will in time create a stable situa- 
tion following successive erosion and slumping of slope material. 

Figure #1 shows a flow net illustrating the character of seepage flow 
through homogeneous and isotropic material in the immediate vicinity of the 
downstream toe of an earth dam, or near the foot of a natural earth slope, 
when immediately underlain by impervious material. The hydraulic gradient 
is equal to the head loss, Ah, represented by any two consecutive equipotential 
curves divided by the distance,41, measured along a flow line, between these 
curves. The direction of seepage flow is normal to the equipotential curves in 
isotropic materials. The seepage therefore is in the direction of the slope 
surface at point “M” where the phreatic line intersects the downstream slope 
while at point “N” the direction of seepage is horizontal. The magnitude of the 
hydraulic gradient at point “M” is equal to sin oc and its magnitude at point 
“N” is tan cc. The seepage force, which represents the effect of fluid friction 
upon the porous material through which it flows, is numerically equal to the 
product of hydraulic gradient by unit weight of water, and has the-same direc- 
tion as the flow. 

Figure # 2 illustrates the forces exerted upon the mineral constituents 
forming two rectangular blocks of soil at the downstream face of the dam, one 
block being located at point “M”, of Figure #1 and the other at point “N”. The 
solid matter is acted upon by its buoyant weight, “W”, and by the seepage 
force exerted by the flowing water, “F”. For equilibrium it is necessary that 
the net force which urges the block of soil down the slope be no greater than 
the available resistance to such motion offered by the soil supporting the block 
on a plane parallel to the face of‘the slope. The buoyant weight of solid matter 
comprising either of these blocks of material is given by the following equa- 
tion: 

S-1 


a W = yw (S-1) (1-n) LD = Ywyz. LD= Yw Sb LD. 


1. Prof., Dept. of Civ. Eng., Univ. of Maine, Orono, Me. 
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wherein e is the void ratio, S, = (S-1) (1-n) is the buoyant specific gravity of 
the solid material as given by the quotient of buoyant unit weight by unit 
weight of water, L and D are the dimensions of the block of earth, and the oth- 
er symbols are defined in Mr. Karpoff’s paper. 

This weight, W, is conveniently decomposed into two components, one nor- 
mal to the slope given by the equation: 


Zz N = W cos a= Yw Sp LD cos a, and 
the other parallel to the slope given by the equation: 
3. T=Wsina= y, S, LD sina. 


The normal component develops friction which may be available in the sup- 
porting soil and hence aids in stabilizing the slope. The tangential component 
tends to move the prism of material down the slope. 

The effect of the seepage force “F” will vary with its direction. In the case 
of point “M” the full magnitude of the seepage force is directed down the 
slope, whereas in the case of point “N” the seepage force has one component 
acting down the slope and a second acting in the opposite direction to the nor- 
mal component of “W”. For these two conditions two different relationships 
determine the maximum angle at which the slope is stable. For the two points, 
“M” and “N”, analysis of the equilibrium of the force vectors of Figure 2 re- 
sults in the relationships derived through equations 4 to 7 inclusive: 


For Point “M” 
F = yy LD sine. 
R =N yw S, LD cosatan®g. 
For Equilibrium it is necessary that 
F+T*R 
or 
Yw LD sina+ yy Sp LD sina® yy, S, LD cosatan¢. 
(1+S,) sina S,- cos atan 


S-1 
Ste 


<Sb_ 
tan as tang 


tan 


For Point “N” 
F = yy LD tana. 
R = (N-F sing) tan$¢= (yy S, LD cosa, LD sinatana) tang. 
For Equilibrium it is necessary that 
F cosa+ TSR 
or 
Yw LD sina=yy S, LD sinaS yy, LD (S,-tan*a) cosatan ¢. 
(1+Sp) sin aS (Sp-tan? a) tan¢cos a. 


-tan? 


i+ Sp tan 


_tan@ <« 1 tan @, or tana 352 
Sp-tan?a 1+Sp 
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Neither of the relationships given in equations 7 agrees with the equation 
22 of Mr. Karpoff’s paper. This is because it is assumed in his analysis that 
the seepage force, while directed parallel to the slope, is given by F = yy LD 
tana. This direction is correct but the magnitude is not for point M, and the 
magnitude is correct but the direction is not for point N. If the values utilized 
in Mr. Karpoff’s numerical example are introduced into the two equations 7, 
one obtains for point “M” a maximum stable slope angle of a= 12° 22 1/2 and 
at point “N” a maximum stable slope angle of a= 11° 51 2/3', as contrasted 
with a slope angle of 12° 14 1/2' by Mr. Karpoff’s method. These values do 
not differ by any great amount and it should be pointed out that as smaller an- 
gles of friction are considered and smaller stable slope inclinations deter- 
mined, the differences in the three methods of computation become smaller 
and smaller. On the other hand, if the angle of friction is increased, then the 
discrepancies between the results of the three methods of determining a will 
also increase. Fig. 3 illustrates graphically the relationship between maxi- 
mum slope angle, a, and angle of internal friction, ¢, in accordance with 
equations 7 and Mr. Karpoff’s equation 22. 

Mr. Karpoff’s equation 25 seems to indicate that the addition of a surface 
layer of highly pervious material will rapidly increase the stability of a given 
slope. It would appear however that if the depth, D, of material which tends to 
move down the slope is very great, the addition of a surface blanket of a given 
thickness, a, would not be as effective in developing stability as it would were 
the dimension, D, very small. Considering the equilibrium of forces at point 
“M” when the face of the slope is surcharged by a free-draining blanket with 
a unit weight of w and a thickness, d, the vertical force acting upon the soil 
prism in the face of the slope is given by equation: 


8. W' = yw Sb LD + w LD. 

and the normal and tangential components of this force by equations: 
9. N' = (yw Sp LD + wLD) cosa. and 

10. T' = (yw Sb LD + wLD) sina. 


The shearing resistance which develops is given by equation: 
11, R' = (yw Sp LD + wLD) cos a@ tan¢. 


and the maximum stable slope angle by equation: 


Sp D + wd 


It is seen that equation 13 involves both the thickness, d, of the blanket and the 
depth, “D”, to which the slope material tends to slide. Figure 4 illustrates, 
for the numerical data used by Mr. Karpoff, the effect on maximum stable 
slope angle, of varying the ratio gf blanket thickness, d, to depth of moving 
slope *D”. It is clear that when D is negligible the blanket is ineffective 


whereas when a becomes very large & approaches @. 

It would appéar that Mr. Karpoff’s equation 25 would indicate the maximum 
slope which would be stable if the depth, “D”, were equal to one foot. The ac- 
tual distance, “D” to which the surface may be unstable is somewhat indeter- 
minate, and it would appear that only by means of a flow net similar to that 
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13. tan a@ = tan d. 


shown in Figure 1 would it be possible, for a given slope, to make a reason- 
able estimate of depth, D, which must be stabilized by the surface blanket. 

In practice the downstream slope of an earthen dam is ordinarily thorough- 
ly protected against the effect of seepage, even if the embankment is con- 
structed of almost homogeneous material. The use of rock toes and horizon- 
tal toe drains to divert the seepage line will, if properly designed, keep the 
seepage line at a substantial distance from the exposed slope. Where the dam 
cross-section is composed of non-homogeneous materials, the downstream 
‘shoulder is ordinarily of more than ample thickness and permeability to en- 
tirely eliminate the need to consider the effect of seepage. 

On the other hand, the construction of slopes in natural ground frequently 
brings about a relocation of the ground water table. Ordinarily no effort is 
made to stabilize such slopes either through reducing the slope angle to a very 
small value or by installing drainage systems. It is frequently possible to di- 
vert the seepage by means of horizontal drains penetrating inward from the 
toe of the slope. In other cases, particularly where soils are fine-grained, 
such drains are not wholly effective unless placed very close together. In 
many instances the expense of installing such drains is excessive. 

The practice of placing a granular blanket on the face of a slope has the ad- 
ditional advantage that it reduces the tendency of the slope to erode under the 
effects of heavy rainfall, and also aids the slope in resisting deterioration 
caused by frost action. It is believed that the treatment of excavated slopes by 
means of such blankets deserves serious consideration in many localities 
where stabilization by other means is not wholly satisfactory. 


JOHN A. FOCHT, Jr.*, A.M. ASCE.—The author has presented a unique 
procedure to approximate the location of the phreatic line for a homogenous, 
isotropic embankment However, a homogenous dam section is unusual in 
present design. Normally, semi-pervious or pervious shells with or without 
drainage blankets are incorporated in the design of the embankment to make 
more effective use of the available borrow materials. In addition, a compact- 
ed fill is probably not isotropic. 

The procedures for analysis of seepage through dams given by Prof. A. 
Casagrande in 1935° enable an engineer to consider both anisotropic and com- 
posite sections. The computation of seepage loss has its greatest potential 
error involved in the determination and selection of the coefficient of permea- 
bility. Therefore, the simplest procedure to estimate the average hydraulic 
gradient through the average area will be satisfactory. A simple flow net can 
be sketched in a very few minutes which will be as accurate as the prediction 
of the embankment permeability. It then appears that the procedure presented 
is not only limited in application but also is more involved than is necessary 
for computation of seepage loss. 

The method for determining the stability of a slope where seepage is 
emerging is correct for cohesionless isotropic material. The factor of safety 
as determined is against the movement of a single particle, but the numerical 
example given is for a clay. No consideration has been given to the cohesive 
strength of the clay. Many dams, levees and dikes built of clay material have 
downstream slopes steeper than 1 vertical to 5 horizontal and have proven to 
be stable. 


2. Associate Eng., Greer & McClelland, Cons. Foundation Engs., Houston, Tex. 
3. “Seepage through Dams”. Arthur Casagrande, Jrnl. of New England Water 
Works Association, Vol. LI, n2, June 1937. 
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The slope in the example for factor of safety of 1.0 would be 1 vertical to 
4.94 horizontal or 12.7 degrees. Common practice for a material with a 
strength of = 23°, C = 0. (which would be a loose, cohesionless silt more 
than a clay) would be to consider a slope of 11.5° (one-half of¢) safe for a 
state of steady seepage. At the same time, the maximum slope which would 
be considered suitable for such material is 23° (¢). The addition of a few feet 
of select granular material on a long slope would not improve the over-all 
stability of the slope against a general failure. It may be conciuded that the 
procedure for checking the stability of the slope surface under steady seepage 
is suitable only for cohesionless materials and not for clays. The mass sta- 
bility of a clay slope must be considered rather than the stability of single 
particles on the slope. 


ABDEL-AZIZ Il. KASHEF*, A.M. ASCE.—Mr. Karpoff presented in his pa- 
per two main parts: the first was the determination of the phreatic line 
through earth dams and the second deals with the slope stability. It is the 
writer’s purpose to discuss only the first of these two parts. 

The seepage of water through dams under the steady state condition is gov- 
erned by the Laplace equation 2. The mathematical solution of this equation 
for such problems, where gravity is the essential driving force, entails some 
difficulty. A very powerful solution, although rather difficult, is the method of 
Hodographs developed by Hamel®. Some mathematical approximations are in 
current use owing to the difficulty involved in exact mathematical solutions*®. 
Among these approximations are Dupuit’s assumptions which are referred to 
in some literature as the Dupuit-Forchheimer Theory. The Pavlovsky’s Theo- 
ry presented by the author is also based on these assumptions. There are 
some objections to applying these assumptions to dam problems as they fail 
to give the true shape of the phreatic line and are not valid for either the en- 
trance or the exit conditions of this line. In water-table well problems where 
a phreatic line exists, it has been concluded? that the failure of the Dupuit- 
Forchheimer Theory to take into cousideration the surface of the line of seep- 
age at the well should alone invalidate its use for determining the shape of the 
free surface. That some engineers are in favor to applying this theory, is be- 
cause of its simplicity and the surprisingly exact value of the flux that can be 
determined by applying it. Experiments °,’ showed that Dupuit’s assumptions 
give the actual distribution of the piezometric heads at the impermeable base 
but do not give the true free surface. Accordingly, Dupuit’s equation is occa- 
sionally referred to as the “base-pressure equation”. Jacob concluded® that 
Forchheimer’s equation for steady-state conditions of water-table wells is 
valid in cases where the ratio of the depth to the lateral extent of the aquifer 
is small and where the range of depth of flow is also not large. He stated that 


4. Lecturer, Ibrahim Univ., Cairo (Egypt); formerly, Research Fellow, Purdue 
Univ., Lafayette, Ind. 

5. Muskat, M., “The Flow of Homogeneons Fluids Through Porous Media,” 
McGraw-Hill, 1937. 

6. Casagrande, A., “Seepage Through Dams,” Journal of New England Water 
Works Association, pp. 131-172, 1937. 

7. Babbitt, H.E. and Caldwell, D.H., “The Free Surface Around and Interfer- 
ence Between, Gravity Wells”, Eng. Expt. Sta., Univ. of Ill., Bull. 374, 1948. 

8. Jacob, C.E., “Flow of Ground Water,” Chapt. V in “Engineering Hydrau- 
lics”, edited by H. Rouse, John Wiley, 1950. 
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these restrictions are, in general no more severe than those that apply, for 
example, to the approximate theory of flexure for beams and concluded: “In 
fact, the engineering theory of beams bears somewhat the same relation to the 
theory of elasticity as the Forchheimer Theory of steady unconfined flow does 
to the more exact theory that may lead to solutions of the Laplace Equation 
with explicit satisfaction of boundary conditions” 

The writer has solved the author’s example by the numerical method—iter- 
ation procedure— which proved to be a very powerful and simple method and it 
does not involve any doubtful assumption. The numerical methods have been 
recently introduced in solving seepage problems. Shaw and Southwell’, in 
1941, solved dam problems by the relaxation method. In 1949, Yang*° solved 
dam and water-table well problems by the same method. The Yang solution 
proved that Dupuit’s assumptions do not give the actual distribution of pieza- 
metric heads at the base. The same conclusion was arrived at by the writer™ 
and others in dealing with water-well problems. The relaxation method was 
also applied by Van Deemter’ in solving a tile drainage problem with a phrea- 
tic line. The relaxation method, though simple, requires good judgment and 
experience. The iteration method, however, is more systematic and far sim- 
pler for an experienced person than the relaxation procedure. In the writer’s 
solution, iteration has been applied in its simplest and lengthy form without 
introducing any expiditing techniques such as “block iteration”**. In spite of 
this, it has been found that the procedure in this way is not laborious. 

A coarse network is laid on the dam cress section, Fig. 1. Any node on 
this network is named as the given numbers of the horizontal and vertical 
lines passing through that node. Thus, for example, the point of intersection 
of the upstream water level with the upstream slope of the dam section is de- 
noted as node 1-18 and so on. A dashed line between two nodes means that any 
of these two nodes will be affected by only one half the value of the other node 
potential. Care should be given to the selection of dashed lines in the upper 
boundary. The initial values of the piezometric heads at the various nodes 
have been roughly assumed and the iteration process has been carried for only 
three traverses. All head values have been multiplied by ten so as to avoid 
writing down the decimals. After this first stage of the coarse network, it has 
been proceeded to a finer net, Figure 2, the size of which is one half the 
coarse net. The values of the new nodes in the finer net are determined from 
the values of the original coarse net nodes as usual.’**!*» “ The upper shape 


9. Shaw, F.S. and Southwell, R.V., “Relaxation Methods Applied to Engineer- 
ing Problems, VII, Problems Relating to Percolation of Fluids Through 
Porous Materials,” Proc. Roy. Soc. London, Series A, vol. 178, pp. 1-17, 
1941. 

10. Yang, S.T., “Seepage Toward a Well Anzlyzed by the Relaxation Method”, 
Doctorate Thesis, Harvard University, 1949. 

11. Kashef, A.I., Touloukian Y.S. and Fadum, R.E., “Numerical Solutions of 
Steady-state and Transient Flow Problems”, Eng. Expt. Sta., Bull. No. 
117, Purdue Univ. U.S.A., July 1952. 

12. Van Deemter, J.J., “Results of Mathematical Approach to Some Flow 
Problems Connected with Drainage and Irrigation”, Applied Scientific Re- 
search, vol. A2, No. 1, pp. 33-53, Netherlands, 1949. 

13. Frocht, M.M., “Photoelasticity”, vol. Il., Chaps. 8 and 9, John Wiley, 
1948. 

14. Grinter, L.E., “Numerical Methods of analysis in Engineering,” McMillan, 
1949. 
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of the finer network is selected according to the deduced free surface of the 
first net. Of course, this first stage could be dispensed with if reasonable ini- 
tial values** were assumed. The phreatic line in the finer net is deter- 
mined"®,"4,5 from the head values at the upper nodes taking into consideration 
that the position head should always equal the piezometric head at this line. 
This condition should also be satisfied for any point at the downstream sur- 
face in the saturated zone above the downstream water level. The iteration 
procedure, for any traverse, has been carried on from upper towards lower 
horizontal lines of the net successively; at each line, the order of successive 
improvements has been taken from left to right. The iteration patterns for 
the various nodes are given in Figures 3 and 4. The iteration method of dif- 
ferences" has been applied by iterating the differences between the head val- 
ues after the first traverse and those of the initial values. In this solution, 
the iterative process terminates after five traverses when it has been found 
that the differences do not exceed * 0.4 ft. The result of the head value at 
node IV-31, Figure 2, is 50.4 ft. while the position head of the point on the 
downstream slope directly above this node is 50.9 ft.; thus point IV-31 may be 
regarded as the discharge point. However, for more accurate determination 
of the phreatic line and the discharge point, one should advance to a more 
finer net which can be done only for part and not for the entire section. 

Table I gives the iteration procedure and results for the upper nodes; 
while Table 2 gives the results of the head values at the other nodes after the 
end of the fifth traverse. The iteration procedure may preferably be made on 
the drawing sheet rather than in a tabulated form. Tables are given in this 
discussion only for convenience. Table 3 includes a comparison between the 
author’s solution using Pavlovsky’s Theory and the iteration solution applied 
in this discussion. The maximum percentage difference, with respect to the 
head difference between the upstream and downstream levels, is about 14.5% 
Although this difference is high, the differences in other examples are of un- 
known magnitudes because this example cannot be regarded as that leading to 
the maximum deviation. 

It is clear that the simplicity and accuracy of the iteration procedure 
should invalidate the application of approximate methods as that presented by 
the author. Moreover, the iteration method determines not only the phreatic 
surface but also the piezometric heads anywhere in the saturated body of the 
dam. It should be noted that iteration procedures do not involve any approxi- 
mations as those of Dupuit’s, and the results become more accurate as the 
nets are selected very fine which may be unwarranted in most of the cases. 
Pavlovsky’s Theory and others cannot be applied to anisotropic soils, in which 
case the iteration method proved"® to be simple and powerful. This shows 
clearly that any theory that involves Dupuit’s assumptions as that presented 
by Mr. Karpoff, can be regarded as merely an approximate method for the de- 
termination of the phreatic line and cannot be applied either for anisotropic 
soils or for the determination of the head distribution in the saturated zone 


below the phreatic line. 


15. Kashef, A.1., “Seepage Through Anisotropic Soils”, Journal of Civil Engi- 
neering, Cairo, April 1954. 
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TABLE 1 


Iteration Solution of the Dam Problem (Fig. 2) 


— Upper Boundary Nedes — 


initial assigned heed values; b : head values 


after first traverse; c = (b - a)3 d, e, f andg: 


4ifferences after end, gre. ang 5™ traverse 


successively; h=(c +d+e+f +g); gj: final 


result of head value = (a = h). 


2817] - 4] - 3] -2]-1 |-2 |-11 | g10 
793] 792] - 1] - 2] -1|-e |- | 785 
II-21 765 | 765 O| -2 |-2 |= 8 | 757 
| 751] 749] - 2|- 7] -1] -2 | -2 |=14 | 737 
| II-23 723) 722| - 1] +6] -1] -2 |-2 © | 723 
II-24] 698] 721] +23/- 1] -2] -3 | -3 |+14 | 712 
III-25 713 | 711] - 2] -13] -7/] -5 | -4 |-31 | 682 
III-26 723| 629] -34] -11] -7] -5 | -4 |-61 | 662 
| 675| 6¢4] -11] - -6] -5 |-3 |-34 | 641 
III-22 645| 637] - 8| - 6| -5 | -4 <3 | -s6 | 619 
ItI-2e | 606] - 3] - 5] ~2] -3 | -2 |-15 | 594 
Iv-3c | 5309] 550] +11] - 3] -2| -1 | 544 
515] - - 2] -i! 504 
| 
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Final Results of the Iteration Solution of the Dam 


Problem (Fig.@) 


"Head Values at the Nodes other than the Upper & 


Lower Boundary ‘Yodes" 


| 

| 

11-16] Iv-15) 247 v-11|878 v-31) 511 745 

11-17] Tv-16| 222 eee 482 722 

II-18|©47]] rv-17}©14 v-13 e658 | v-33] 463 

rv-19] 728 v-14]247 | v-34| 453 v1-23] 601 

III-14| rv-19] 79¢ v-15] 233 vr- 4] 900 vr-241677 

III-15]£6°]| Iv-20| 760 VI- 5] 897 659 

rv-21| 742 |] v-17] 7¢¢ 895 |] vt-26] 643 

III-17| Iv-22| 726 784 | vI- 7/291 623 j 

rII-12|218 Iv-z23] 710 v-1¢] 768 | vI- sso see 

III-1¢] tv-24] cs2 750 | VI- vr-2e] 572 

111-20| 772] 1v-25] 675 v-21| 732 vr-10] 876 || vz-30| 547 q 

111-211 751] rv-26] 657 v-2e] 716 vr-11] see |} vz-31] 527 

111-22|732] 1v-27| v-23| 702 vr-12| aso 492 

III-23]719]] Iv-28 v-24| 684 |] vr-13] 849 vr-33] 473 
1v-29] v-25| ces vi-14| eae v1-34] 462 

1v-10]205]] Iv-32| 455 v-26|650 v1-15] 823 vr-35| 455 

v-7 v-27| 622 vi-16| ece v1-36] 453 | 

Tv-12|ge2 v-8 vr-17| 791 vr-37] 452 

Iv-13|°73]] v-2e|576 || vr-18| 776 
Iv-14]e63]] v-10]eg2]] v-30] 545 |] vr-19] 762 

573-12 

j 

al 


TABLE 3 


Comparison between Pavlovsky’s Theory and 
Iteration 3olutions for the Dam 
Problem ( Fig. 2 ) 


=7.0 ft. (Pav.) 162 ft. (Pav.) 
5.4 ft. (Iter.) 165 ft. (Iter.) 


Piezometric Head Difference 


Pav. (in feet)}| Iter. (in feet) | (in feet) 


PsSa? a , x and S have the same significance as 


in the original paper. 


> 
Bo 
72.0 78.5 6.5 
15 70 75.5 5.1 
30 68 8 73.7 4.2 
; 45 67.1 72.3 5.2 
60 65.4 71.2 5.8 
75 63.5 68.2 4.7 
20 61.7 66.2 4.5 
120 57.9 61.¢ 4.0 
| 135 55.8 59.4 3.6 
150 53.7 54.4 0.7 
| 
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WILLIAM H. GOTOLSKI,"° J.M. ASCE. - The recognized importance of de- 
termining the location of the phreatic surface in earthen hydraulic structures 
makes it necessary that the formulas and procedure covering the location of 
such surfaces be of ready application. It is realized that the exact determina- 
tion of the position of the phreatic surface is not necessary; as a matter of 
fact, present knowledge as to the nature of the action of capillary water in this 
area does not yet permit making precise determinations. 

The method presented by Mr. Karpoff for locating the phreatic surface is 
interesting and simple to use. The writer wishes to supplement Mr. Karpoff’s 
method with a point-by-point comparison between the author’s method and an- 
other procedure, after Robert Dachler,'’ wl.ich seems to have similar advan- 
tages of speed and simplicity. 

It will be noted that this discussion deals only with embankments on imper- 
vious foundations because from the hydraulic standpoint an embankment on an 
impervious foundation results in the highest position of the phreatic line which 
in turn tends to decrease the stability of the slope. It is also pointed out that 
the tailwater depth will be considered negligible here since the tailwater depth 
in practice is usually zero. 

The development of formulas for seepage flow through an earthen dam and 
for locating the phreatic surface in Mr. Dachler’s method and the author’s 
method are similar in certain respects. In both cases, for purposes of sim- 
plification, an embankment of homogeneous material and of unit width, is di- 
vided into three parts. The three characteristic parts of Mr. Karpoff’s meth- 
od may be found in the author’s Figure 1. Mr. Dachler’s cross-section is to 
be found in the writer’s Figure 1. These separate flow systems, into which 
the dam cross-section is separated by different approximations, represents 
one hydraulic unit, and through each portion the same amount of seepage, q, 
occurs. It is also noted that the fluid heads should be continuous in passing 
from one to the other portion. Both methods agree on this requirement. 

In making a point-by-point comparison of the two methods, let us proceed 
from the upstream part or inflow part of the dam. The author’s method con- 
siders the upstream part as being that portion of the dam between the upstream 
slope and Section W-W (see author’s Figure 1) and the streamline are as- 
sumed therein to be straight and horizontal. The Section W-W was chosen as 
a boundary of this portion because of geometric simplicity and laboratory ex- 
periments which indicated that the transition of the flow line from this up- 
stream part to the central part of the dam occurred in the vicinity of this sec- 
tion.!® The equation for this portion of the dam may be found in the author’s 
paper. 

Mr. Dachler assumed the inflow section of the system to be bounded by 
ABCD in the writer’s Figure 1. The flow through this region is computed by 
means of a conjugate-function transformation suitable for a wedge of angle A, 
bounded by streamline AD and equipotential line AB. BC is a streamline leav- 
ing B, with no attempt to correct for the fact it is not a free surface stream- 
line. 

If A = 277/n, n being an integer, the complex-variable transformation that 
is suitable may be expressed as 


16. Instructor of Civ. Eng., The Pennsylvania State Univ., State College, Pa. 

17, “Grundwasserstroemung,” by Robert Dachler, Julius Springer, Vienna, 
1936. 

18. U.S. Bureau of Reclamation Translation, Tech. Memo. No. 383. June 1934. 
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H+i W = (x + iy)", (1) 


where the velocity potential has been replaced by the fluid head H. This sub- 
stitution is for convenience only. 

The location of the potential CD can be made definite by requiring the loss 
of energy by friction in passing through ABCD be equal to that loss in passing 
through ABE. This loss of energy experienced by the fluid in passing between 
two streamlines separated by AWacross the drop AH in the fluid head is 


dw = Yg AHAY, (2) 


where y is the fluid density and g is the acceleration of gravity. The result 
of this condition becomes equivalent to the requirement that BCF and FDE 
have the same area in the (H, W) plane.’® By use of conformal mapping the 
resultant flux of ABCD is found to be expressed as a function of the angle A 
and is called the form factor: 

f, = = 1.12 + 1.93tana (3) 


It is to be noted that a in Equation (3) refers to the potential drop only over the 
segment 8C of the free surface. 

The fundamental difference between these methods of attacking the problem 
of the flow in the upstream portion of the dam lies in that fact that the author’s 
method assumes the streamlines to be horizontal and straight between the up- 
stream slope and Section W-W, whereas Dachler uses a transformation from 
an area bounded by flow lines and potential lines to a triangular area for pur- 
poses of simplification. Both methods have good points in their favor for such 
assumptions. It is not the purpose of this discussion to argue the case for ei- 
ther method. 

The flow to the right of BE is further resolved into separate parts, the main 
body of the dam and the downstream part. The usual case for the downstream 
part of the dam is considered here, i.e., zero tailwater. 

In locating the emergence point consider the branch M,GC (writer’s Fig- 
ure 2). It is assumed that the flow in this branch is horizontal. Applying Du- 
puit’s fundamental equation, and constructing the potential line M,C as a para- 
bola, this permits the segment length L , to be expressed as a f(G, yo, y). 
Consequently 


Le-y(l -—_) tanG + ycotG (4) 
2Yo 


The total discharge q may then be determined by integration and, reduced 
through a series formation, may be expressed as 
sin2G _q _ sin2G 

In determining the point of emergence in the author’s method, it is assumed 
that the flow is horizontal in the downstream part (the portion between the 
downstream slope and Section W,, see author’s Figure 9). It may be noted 
that both methods assume the flow to be horizontal in this portion. The differ- 
ence in the location of the point of emergence appears to vary due to the shape 
of the potential line which passes through the point of emergence, M,. Pavlov- 
sky assumes the line to be straight and Dachler uses a parabolic curve. It 


19. “The Flow of Homogeneous Fluids Through Porous Media,” by Morris 
Muskat, McGraw-Hill Book Co., In., New York, N.Y., 1937. p. 340. 
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appears that the assumption of a parabolic potential line is much closer to the 
actual flow case as the phreatic line itself is a parabola. It follows then, from 
the equipotentiai theory, that the potential lines should be perpendicular to the 
flow lines. Consequently a parabolic potential line would fill this requirement 
more closely than would a straight line. 

Dachler considered the flow in the main body of the dam to be essentially 
linear in character, the flow being related to the drop in fluid head across it, 
H,, by 


_ _ Hp 
4° kH, (6) 

where Hyp is the mean height of the phreatic line along CM,, and L is the dis- 

tance EN (see writer’s Figure 1). This form factor is used to relate the main 

body of the dam to the inflow and downstream portions of the dam. 

The extent of the central part or main body of the dam differs in the au- 
thor’s and writer’s methods by the area bounded by Section W-W (author’s 
Figure 1) and BE (writer’s Figure 1). This is due to the initial assumption of 
the make up of the inflow or upstream portion of the dam. In both Mr. Kar- 
poff’s method and Mr. Dachler’s,” by applying Darcy’s law for the flow of wa- 
ter through homogeneous material and Dupuit’s relation, the general equation 
for the phreatic line in the central part or main body of the dam may be de- 
rived. According to Dachler the expression is 


where q = seepage flow per unit length of dam 
k*? = coefficient of permeability 
L = horizontal distance from the point of emergence, M,, to the desired 
point on the phreatic line 
y, = vertical ordinate of point M,, above the dam base 
y, = desired vertical ordinate above dam base. 
Consequently the procedure to be used is as follows: First, angle G be- 
ing known, assume a discharge q. By using Equation (5), y, is computed. 
Thus the point of emergence is temporarily located. 

Then using the assumed q as above and solving Equation (3), a is deter- 
mined, Since a + H, + y, = H, H, may be found. Knowing H,, compute Hm. 
Solve Equation (6) for q; lack of. agreement with the assumed value of q re- 
quires a repetition of this process until agreement is reached. 

With q, a, and y, knowa, it is now possible to determine the position of the 
phreatic line in the main portion of the dam. Equation (7) is used for that 
purpose, and in the following form: 


2qL 
yy? (8) 


It is pointed out during this discussion that there are several pertinent dif- 
ferences in these methods of locating the phreatic surface in adam. The 


20. “Grundwasserstroemung,” by Robert Dachler, Julius Springer, Vienna, 
1936, p. 90. 

21. It is assumed for convenience here and throughout this discussion that the 

effective permeability k has the value of unity. 
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differences in the profiles indicate that the author’s method results in large 
deviations from other methods. For further actual comparison and for illus- 
trative purposes the following example is used. 


Illustrative Example 
Given: Height of dam - 100.0' 
Freeboard 
Headwater level - 
Width of crest 
Upstream slope 
Downstream slope - 
Find: 
(1) the point of emergence of the phreatic line on the downstream slope 
(2) the maximum ordinate at the intersection of the phreatic line with a. 


To obtain a solution to this problem reference is made to the procedure 
outlined in the previous paragraphs. A q is assumed; determinations of a, H, 
and y, are made. A check is then made on q. This method results in the re- 
quired values of q, a, H, and y,. 

A q of 11.5 was assumed and found to check out by use of the procedure 
outlined above. Thus y, was found to be 38.4 ft. Solving Equation (3) for a re- 
sulted in 6.5 ft. H, is equal to 90 minus 44.9 or 45.1 ft. Solving for Hm re- 
sulted in a check of the assumed q. Consequently the point of emergence oc- 
curs on the downstream slope at a point which is 38.4 ft above the dam base. 
Starting at this point and using Equation (8) the equation of the curve is 


2 _2X11.5L 
“ k 


where L is the distance from the point of emergence to the desired point on 
the phreatic line. The plotted data are shown in the writer’s Figure 3. 

Table I is a comparison showing the ordinates of the phreatic line obtained 
by various methods. It is noted that the author’s method results in a lower 
position of the phreatic line than do any of the other methods plotted. It would 
seem then that the author’s method should be questioned as to its use for mak- 
ing stability analyses. The writer is of the opinion that the higher the position 
of the phreatic line the less stable is the slope. It is thus concluded that the 
method giving the highest position of the phreatic line should be used when 
slopes are to be analyzed for stability. 

To complete the phreatic line the reversed curve BC (see writer’s Figure 1) 
is drawn perpendicular to the upstream slope at the entrance point and tangent 
to Dupuit’s curve at the point where Dupuit’s curve intersects a. This follows 
from the equipotential theory and the continuity of the flow. 

In conclusion, the writer wishes to emphasize his belief that the assump- 
tions of Mr. Dachler approach the actual conditions more closely than do those 
of the author’s method. However, Mr. Karpoff has made a valuable contribu- 
tion in presenting his paper, especially the mathematical approach which may 
be used by field engineers to make quick decisions for determining the stabil- 
ity of an earthen dam. 


+ (38.4) (9) 
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Table I—The values of x and y. 


Coordinates 


Note: The origin is taken as the toe of the downstream slope. 
*y, or height of point of emergence above the dam base. 


y in feet 
X, ft DACHLER | L. CASAGRANDE | A. CASAGRANDE |} PAVLOVSKY 
85 ~ - - 27.85* | 
102.5 34* 33 
109.5 _ 37 36.5 35 
115 38.4* 38 38 36 
160 51 51 51 46.5 
200 60 60 60 54.5 | 
240 66.5 67.5 66.5 60.5 
280 73.5 714 71.5 65.5 
f 320 80 80 76.5 70.5 
330 81.5 82.5 77.5 72.5 
340 84 83.5 80.5 76.5 
350 88 88 88 88 
| 
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SAND DRAIN APPLICATIONS BY THE 
PORT OF NEW YORK AUTHORITY 


by John M. Kyle and Martin S. Kapp. 
(Proc. Sep. 456) 


K.B. HIRASHIMA.*—The authors have presented valuable data on field ex- 
perience with sand drains. It is evident from reading the paper that all of the 
work described was carefully planned and expertly executed. 

The Territorial Highway Department has used the sand drain technique on 
two highway projects and the results obtained were presented before the last 
meeting on the Highway Research Board.” On a third project, sand drains 
were considered but were later omitted due to a change in alignment. 

The Department’s previous experience with highways constructed over 
swamp areas had been unsatisfactory. The settlements continued for years, 
since about 1937 for one project, with no end in sight. This has made it neces- 
sary for the maintenance forces to go back every few years, re-level and 
bring the road surface up to grade. The primary objective in using sand 
drains on the above two projects was, therefore, to accelerate the ultimate set- 
tlement, causing the latter to take place within a reasonable time interval. In 
this respect, the use of sand drains proved successful. Both projects were 
completed rather quickly, before the full theoretical settlement had taken 
place. But the additional settlements that have taken place since completion 
of the projects have been relatively small. 

On both of the Department’s sand-drain projects, the actual settlements 
with time were carefully recorded and the writer has found it interesting to 
compare the results with those of the authors. 

It is the Department’s experience that there is a distinct increase in rate 
of settlement immediately upon installation of sand drains. But there is a still 
greater rate of increase as the embankment and surcharge loads are applied. 

As the embankment load increases, the pore-water pressure also increases. 
Usually the permeability of the soil is low so that the pore-water ordinarily 
cannot drain readily. Under such conditions, if the loading is increased too 
rapidly, a dangerously high internal hydraulic pressure will be built up, which 
can result in mud flows. The function of sand drains is to provide additional 
drainage channels so as to more quickly dissipate the high internal pore-water 
pressures. Even then, on contract jobs, the tendency is to apply the load too 
rapidly rather than too slowly. As a consequence, part of the settlement dur- 
ing the loading period in such cases is due, not to drainage of the pore-water, 
but to plastic displacement. However, after a short time interval, the plastic 
displacement will cease, if the overloading was moderate and practically all 
the subsequent settlement will be due to drainage. If the rate of loading is 


* Materials Eng., Territorial Highway Dep., Honolulu, Hawaii 
22. “Hawaii’s Experience with Vertical Sand Drains” by K.B. Hirashima, a pa- 
per presented at the 33rd Annual Meeting of the Highway Research Board. 
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very rapid, complete failure by lateral displacement of the undersoil, includ- 
ing the sand drains, is possible. This is a fault not of sand drains as such but 
of a faulty technique. 

In all of the settlement curves presented by the authors, the actual settle- 
ments appear to be in good agreement with those calculated. Exact agreement 
is hardly to be expected in the case of soils, which as is well-known, can vary 
widely in physical properties from point to point. 

Figure A shows a comparison of the calculated with the actual observed 
field settlements taken from one of the Department’s projects. Prior to the 
installation of sand drains, a blanket of fill material was placed over the 
swamp to serve as a working table for the construction equipment. The time 
when construction work was started on the working tabie is taken as zero on 
the time axis. The observed settlements necessarily reflect the fact that the 
full load was not applied suddenly but over a period of time. Accordingly, for 
the calculated curve, the load was assumed to be applied in increments and 
not instantaneously. As far as long-term effects are concerned, say settle- 
ments at time exceeding six months or so after application of full load, the re- 
sults are practically the same whether the loading is assumed to be instanta- 
neous or gradual. However, for near-term effects, that is the settlements dur- 
ing the first few weeks after installation of sand drains, the difference is dis- 
cernable. If the calculated curve had been computed on the basis of instanta- 
neous loading, the apparent agreement between the calculated and observed 
curves would have been closer in this case, although Figure A is on too small 
a scale to show this. In connection with the Newark Airport project, the au- 
thors mentioned slides caused by too rapid loading. This was also true to 
some extent in the case illustrated by Figure A. From the curves, the plastic 
displacement caused by this rapid loading appears to have been about 0.5 ft. 
That is, if this plastic displacement had not occurred, the agreement between 
observed and calculated settlements would have been closer. 

According to the data presented by the authors, the ratio of the horizontal 
to vertical permeability was low for all the projects covered, not exceeding | 
two. Presumably, the values were determined by laboratory tests. The De- 
partment’s experience is that such values are more easily and reliably deter- 
mined by means of rather simple field tests. For the vertical permeability, 
the “tube” method? was used. A good approximation of the horizontal perme- 
ability is given by the so-called “auger” hole method.™%»* Both methods are 
described in “Soil Engineering” by Spangler. However, for the auger hole me- 
thod, the constants given in Spangler’s book appear to the writer to give values 
that are too high. , 

The value of the paper would have been enhanced if cost data on sand drains 
had been included. On the Department’s projects, the sand drains were 18" 
in diameter and the costs ranged from $1.05 to slightly over $1.26 per linear 
foot.?? 


23. “A Field Method for Measuring the Permeability of Soil Below A Water 
Table” by Richard K. Frevert and Don Kirkham; Proc. Highway Research 
Board, 1948, Vol. 28, pp. 433-442. 

24. “Theory of Seepage Into Auger Holes” by Don Kirkham and C.M.H. Van 
Bavel; Proc. Soil Science Society of America, 1948, pp. 75-82. 

25. “Field Measurement of Soil Permeability Using Auger Holes” by C.M.H. 
Van Bavel and Don Kirkham; Proc. Soil Science Sciety of America, pp. 90- 
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PROCEEDINGS-SEPARATES 


The technical papers published in the past year are presented below. Technical-division sponsorship is indicated by an 
abbreviation at the end of each Separate Number, the symbols referring to: Air Transport (AT), City Planning (CP), Con- 
struction (CO), Engineering Mechanics (EM), Highway (HW), Hydraulics (HY), Irrigation and Drainage (IR), Power (PO), 
Sanitary Engineering (SA), Soil Mechanics and Foundations (SM), Structural (ST), Surveying and Mapping (SU), and Water- 
ways (WW) divisions. For titles and order coupons, refer to the appropriate issue of “Civil Engineering” or write for a 
cumulative price list. 


VOLUME 79 (1953) 


DECEMBER: 359(AT), 360(SM), 361(HY), 362(HY), 363(SM), 364(HY), 365(HY), 366(HY), 367(SU)°, 368(WW)°, 369(IR), 
370(AT)©, 371(SM)©, 372(CO)©, 373(ST)©, 374(EM)°, 375(EM), 376(EM), 377(SA)°, 378(PO)°. 


VOLUME 80 (1954) 


JANUARY: 379(SM)°, 380(HY), 381(HY), 382(HY), 383(HY), 384(HY)°, 385(SM), 386(SM), 387(EM), 388(SA), 389(SU)°, 390(HY), 
391 (IR), 392(SA), 393(SU), 394(AT), 395(SA)°, 396(EM)°, 397(ST)°. 


FEBRUARY: 398(IR)4, 399(SA)4, 400(CO)4, 401(SM)°, 402(AT)4, 403(AT)4, 404(IR) 4 405(P0)4, 406(AT)4, 407(su)4, 408(sU)4, 
409(ww)4, 410(AT)4, 411(SA)9, 412(PO)4, 413(HY)4. 


MARCH: 414(ww)4, 415(su)4, 416(sm)°, 417(SM)4, 418(AT)4, 419(SA)4, 420(SA)4, 421(AT)4, 422(SA)4, 423(CP)4, 424(AT)4, 
425(SM)4, 426(IR)4, 427(ww)4. 


APRIL: 428(HY)°, 429(EM)°, 430(ST), 431(HY), 432(HY), 433(HY), 434(ST). 
MAY: 435(SM), 436(CP)°, 437(HY)°, 438(HY), 439(HY), 440(ST), 441(ST), 442(SA), 443(SA). 


JUNE: 444(SM)®, 445(SM)®, 446(ST)®, 447(ST)®, 448(ST)®, 449(ST)®, 450(ST)®, 451(ST)®, 452(SA)®, 453(SA)®, 454(SA)®, 
455(SA)®, 456(SM)®. 


JULY: 457(AT), 458(AT), 459(AT)°, 460(IR), 461(IR), 462(IR), 463(IR)°, 464(PO), 465(PO)°. 


AUGUST: 466(HY), 467(HY), 468(ST), 469(ST), 470(ST), 471(SA), 472(SA), 473(SA), 474(SA), 475(SM), 476(SM), 477(SM), 
478(SM)°, 479(HY)°, 480(ST)°, 481(SA)°, 482(HY), 483(HY). 


SEPTEMBER: 484(ST), 485(ST), 486(ST), 487(CP)° 488(ST)°, 489(HY), 490(HY), 491(HY)°, 492(SA), 493(SA), 494(SA), 
495(SA), 496(SA), 497(SA), 498(SA), 499(HW), 560(HW), 501(HW)°, 502(WW), 503(WW), 504(WW)°, 505(CO), 506(CO)°, 
507(CP), 508(CP), 509(CP), 510(CP), 511(CP). 


OCTOBER: 512(SM), 513(SM), 514(SM), 515(SM), 516(SM), 517(PO), 518(SM)°, 519(IR), 520(IR), 521(IR), 522(IR)°, 523(AT)°, 
524(SU), 525(SU)°, 526(EM), 527(EM), 528(EM), 529(EM), 530(EM)°, 531(EM), 532(EM)°, 533(PO). 


NOVEMBER: 534(HY), 535(HY), 536(HY), 537(HY), 538(HY)°, 539(ST), 540(ST), 541(ST), 542(ST), 543(ST), 544(ST), 545(SA), 
546(SA), 547(SA), 548(SM), 549(SM), 550(SM), 551(SM), 552(SA), 553(SM)°, 554(SA), 555(SA), 556(SA), 557(SA). 


DECEMBER: 558(ST), 559(ST), 560(ST), 561(ST), 562(ST), 563(ST)°, 564(HY), 565(HY), 566(HY), 567(HY), 568(HY)°, 569(SM), 
* §70(SM), 571(SM), 572(SM)°, 573(SM)°, 574(SU), 575(SU), 576(SU), 577(SU), 578(HY), 579(ST), 580(SU), 581(SU), 582(Index). 


c. Discussion of several papers, grouped by Divisions. 
d. Presented at the Atlanta (Ga.) Convention of the ‘ciety in February, 1954. 
e. Presented at the Atlantic City (N.J.) Convention in June, 1954. 


| 

: 

| 

| 

| 

7 


AMERICAN SOCIETY OF CIVIL ENGINEERS 


OFFICERS FOR 1955 


PRESIDENT 
WILLIAM ROY GLIDDEN 


VICE-PRESIDENTS 


Term expires October, 1955: Term expires October, 1956: 
ENOCH R. NEFDLES FRANK L. WEAVER 
MASON G. LOCKWOOD LOUIS R. HOWSON 


DIRECTORS 


Term expires October, 1955: Term expires October, 1956: Term expires October, 1957: 
CHARLES B. MOLINEAUX WILLIAM S. LaLONDE, JR. JEWELL M. GARRELTS 
re ge OLIVER W. HARTWELL FREDERICK H. PAULSON 
CARL G. PAULSEN THOMAS C. SHEDD GEORGE S. RICHARDSON 
LLOYD D. KNAPP SAMUEL B. MORRIS DON M. CORBETT 
GLENN W. HOLCOMB ERNEST W. CARLTON GRAHAM P. WILLOUGHBY 
FRANCIS M. DAWSON RAYMOND F. DAWSON LAWRENCE A. ELSENER 


PAST-PRESIDENTS 


Members of the Board 


WALTER L. HUBER DANIEL V. TERRELL 


EXECUTIVE SECRETARY TREASURER 
WILLIAM N. CAREY CHARLES E. TROUT 


ASSISTANT SECRETARY ASSISTANT TREASURER 
E. L. CHANDLER CARLTON S. PROCTOR 


PROCEEDINGS OF THE SOCIETY 


HAROLD T. LARSEN 
Manager of Technical Publications 


DEFOREST A. MATTESON, JR. PAUL A. PARISI 
Editor of Technical Publications Assoc. Editor of Technical Publications 


COMMITTEE ON PUBLICATIONS 
SAMUEL B. MORRIS, Chairman 
JEWELL M. GARRELTS, Vice-Chairman 
GLENN W. HOLCOMB OLIVER W. HARTWELL 
ERNEST W. CARLTON DON M. CORBETT 


4 
7 
- 
q 
| 
| 
a 
+ 
j ? 
q 


